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We construct monopole-antimonopole chain and vortex solutions in Yang-Mills-Higgs theory cou- 
pled to Einstein gravity. The solutions are static, axially symmetric and asymptotically flat. They 
are characterized by two integers (m, n) where m is related to the polar angle and n to the azimuthal 
angle. Solutions with n = 1 and n — 2 correspond to chains of m monopoles and antimonopoles. 
Here the Higgs field vanishes at m isolated points along the symmetry axis. Larger values of n 
give rise to vortex solutions, where the Higgs field vanishes on one or more rings, centered around 
the symmetry axis. When gravity is coupled to the flat space solutions, a branch of gravitating 
monopole-antimonopole chain or vortex solutions arises, and merges at a maximal value of the cou- 
pling constant with a second branch of solutions. This upper branch has no flat space limit. Instead 
in the limit of vanishing coupling constant it either connects to a Bartnik-McKinnon or generalized 
Bartnik-McKinnon solution, or, for m > 4, n > 4, it connects to a new Einstein- Yang-Mills solution. 
In this latter case further branches of solutions appear. For small values of the coupling constant 
on the upper branches, the solutions correspond to composite systems, consisting of a scaled inner 
Einstein- Yang-Mills solution and an outer Yang-Mills-Higgs solution. 

PACS numbers: 04.20.Jb, 04.40.Nr 

I. INTRODUCTION 

The non-trivial vacuum structure of SU(2) Yang-Mills-Higgs (YMH) theory allows for the existence of regular 
non-perturbative finite energy solutions, such as monopoles, multimonopoles and monopole-antimonopole systems. 
While spherically symmetric monopoles carry unit topological charge 0, , monopoles with charge n > 1 are axially 
symmetric @, 0,M El ° r possess no rotational symmetry at all 0, U • 

In the Bogomol'nyi-Prasad-Sommerfield (BPS) limit of vanishing Higgs potential, the monopole and multimonopole 
solutions satisfy a set of first order equations, the Bogomol'nyi equations 0. The spherically symmetric and axially 
symmetric BPS (multi)monopole solutions are known analytically 0, 0- In these solutions all nodes of the Higgs 
field are superimposed at a single point. In BPS multimonopole solutions with only discrete symmetries, recently 
constructed numerically, the nodes of the Higgs field can be located at several isolated points The energy of BPS 
solutions satisfies exactly the lower energy bound given by the topological charge. In particular, in the BPS limit 
the repulsive and attractive forces between monopoles exactly compensate, thus BPS monopoles experience no net 
interaction [To| . 

The configuration space of YMH theory consists of topological sectors characterized by the topological charge of the 
Higgs field. As shown by Taubes each topological sector contains besides the (multi)monopole solutions further 
regular, finite energy solutions, which do not satisfy the first order Bogomol'nyi equations, but only the set of second 
order field equations. These solutions form saddlepoints of the energy functional [ilj. The simplest such solution, the 
monopole-antimonopole pair solution [l2l 1 1 3l | , is topologically trivial. It possesses axial symmetry, and the two nodes 
of its Higgs field are located symmetrically on the positive and negative z-axis. 

Recently we have constructed axially symmetric saddlepoint solutions, where the Higgs field vanishes at m > 2 
isolated points on the symmetry axis. These represent chains of single monopoles and antimonopoles located on the 
symmetry axis in alternating order For an equal number of monopoles and antimonopoles, i.e. for even m, the 
chains reside in the topologically trivial sector. When the number of monopoles exceeds the number of antimonopoles 
by one, i.e. for odd m, the chains reside in the sector with topological charge one. While such chains can also be 
formed from doubly charged monopoles and antimonopoles (n = 2), chains with higher charged monopoles (n > 2) 
cannot be formed. Instead completely different solutions appear, where the Higgs field vanishes on one or more rings 
centered around the symmetry axis |l4j . We therefore refer to these solutions as vortex rings. 

When gravity is coupled to YMH theory, this has significant effect on the monopole, multimonopole and monopole- 
antimonopole pair solutions. When the coupling constant is increased from zero |l5l Il6j . a branch of gravitating 
monopole solutions emerges smoothly from the flat space 't Hooft-Polyakov monopole solution. The coupling constant 
a, entering the Einstein- Yang-Mills-Higgs (EYMH) equations, is proportional to the Higgs vacuum expectation value 
77 and the square root of the gravitational constant G and to the square of the Higgs vacuum expectation value r\. With 
increasing coupling constant a the mass of the gravitating monopole solutions decreases. The branch of gravitating 
monopole solutions extends up to a maximal value a max , beyond which gravity becomes too strong for regular 



2 



monopole solutions to persist |15|, For vanishing Higgs self-coupling constant, this first gravitating monopole 
branch merges with a second branch at a max > which extends slightly backwards, until at a critical value a cr of the 
coupling constant a degenerate horizon develops ^3] . The exterior space time of the solution then corresponds to the 
one of an extremal Reissner-Nordstrom (RN) black hole with unit magnetic charge At a cr this second branch 
of gravitating monopole solutions thus bifurcates with the branch of extremal Reissner-Nordstrom (RN) solutions. 
Remarkably, the additional attraction in the YMH system due to the coupling to gravity also allows for bound 
monopoles, not present in flat space [lq . 

For the monopole-antimonopole pai r solution, on the other hand, we observe a different coupling constant depen- 
dence for the gravitating solutions |l8j . Again a branch of gravitating monopole-antimonopole pair solutions emerges 
smoothly from the flat space solution, and merges at a maximal value of the coupling constant a max with a second 
branch of gravitating monopole-antimonopole pair solutions. This second branch, however, extends all the way back 
to vanishing coupling constant. Along the first branch the mass of the solutions decreases with increasing a, since with 
increasing gravitational strength the attraction in the system increases. Along the second branch the mass increases 
strongly with decreasing coupling constant and diverges in the limit of vanishing coupling constant. This upper 
branch may be considered as being obtained by decreasing the Higgs expectation value. Along the upper branch the 
solutions shrink to zero size in the limit a — > 0. Scaling the coordinates, the mass and the Higgs field by a shows, that 
the scaled solutions approach the lowest mass Bartnik-McKinnon (BM) solution of Einstein- Yang-Mills (EYM) 
theory in the limit a — > 0. 

Here we consider the effect of gravity on the monopole-antimonopole chains and vortex rings [hH ] . Characterizing 
these solutions by the integers m and n, related to the polar angle 9 and the azimuthal angle (p, respectively, we 
first consider chains, where n < 2, and then vortex rings, where n > 3. When m is even, these solutions reside in 
the topologically trivial sector, thus we expect a similar dependence on the coupling constant as for the monopole- 
antimonopole pair. When m is odd, the solutions reside in the sector with charge n, where the corresponding 
multimonopoles reside. Therefore at first sight a connection with RN solutions appears to be possible, similar as 
observed for gravitating (multi)monopoles. We find, however, that all these unstable gravitating solutions, chains and 
vortex rings alike, show the same general coupling constant dependence as observed for the monopole-antimonopole 
pair. Two branches of solutions arise, a lower branch connected to the flat space solution and an upper branch 
connected to an EYM solution. For m < 3 these EYM solutions correspond to the spherically symmetric BM 
solutions (n= 1) or their axially symmetric generalizations (n > 1) [Hi l20j|. For m > 4 new EYM solutions arise in 
the limit |2l|, implying further branches of gravitating EYMH solutions. 

In this paper we present gravitating monopole-antimonopole chains and vortex rings for vanishing Higgs self- 
coupling. In section II we present the action, the axially symmetric Ansatz and the boundary conditions. In section 
III we discuss the properties of these solutions with particular emphasis on the limit a — > 0, and we present our 
conclusions in section IV. 



II. EINSTEIN- YANG-MILLS-HIGGS SOLUTIONS 



A. Einstein- Yang-Mills-Higgs action 

We consider static axially symmetric monopole-antimonopole chains and vortex rings in SU(2) Einstein- Yang-Mills- 
Higgs theory with action 

with curvature scalar R, su(2) field strength tensor 

Ffi V = d^Av - dyA^ + ie[A„, A v ] , (2) 

gauge potential A^ = A^T a /2, and covariant derivative of the Higgs field $ = $ a r a in the adjoint representation 

D^ = d^ + ie[A fl ,$] . (3) 

Here G and e denote the gravitational and gauge coupling constants, respectively, 77 denotes the vacuum expectation 
value of the Higgs field, and A the strength of the Higgs self-coupling. 
Under SU(2) gauge transformations U, the gauge potentials transform as 

A; = UApUi + -(d^U)^ , (4) 
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and the Higgs field transforms as 

= . (5) 

The nonzero vacuum expectation value of the Higgs field breaks the non-Abelian SU(2) gauge symmetry to the 
Abelian U(l) symmetry. The particle spectrum of the theory then consists of a massless photon, two massive vector 
bosons of mass M v = er/, and a massive scalar field M s = \[2\r\. In the BPS limit the scalar field also becomes 
massless, since A = 0, i.e. the Higgs potential vanishes. 

Variation of the action with respect to the metric g^ u leads to the Einstein equations 

G>„ = i? p „ — -g^R = 8irGTfi„ (6) 

with stress-energy tensor 

rr, T o 0L M 

= 2 Tr (F^F^ - \ 9llv F a pF^) 

+ Tr {\D^D^ - ± 9lxv D a <f>D<*$) - ^„Tr(<I> 2 - if) 2 . (7) 
Variation with respect to the gauge field and the Higgs field $ leads to the matter field equations, 

1 :D^J—gF^) - iie[$, £>"$] - , (8) 



-g * ' 4 
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respectively. 



DpiV^gD**) + \(& - r? 2 )$ = , (9) 



B. Static axially symmetric Ansatz 



To obtain gravitating static axially symmetric solutions, we employ isotropic coordinates |Tft 12(1 l2lj. In terms of 
the spherical coordinates r, 9 and ip the isotropic metric reads 

, o m , mr 2 ,„, lr 2 sin 2 9 , 9 . 
ds 2 = -fdt 2 + jdr 2 + —d0 2 + - dip 2 , (10) 

where the metric functions /, m and / are functions of the coordinates r and 9, only. The z-axis (9 = 0, 7r) represents 
the symmetry axis. Regularity on this axis requires \22\ 

m\$=o,-K = l\e=o,ir • (11) 
We parametrize the gauge potential and the Higgs field by the Ansatz 01 



Kl . , r-> Ty \ ir\ T V 



A^dx" = —dr + (1 - K 2 )d9 



(n) 



((n,m) (n,m) \ 

+ ( l - K ^\^\ d( P ' ( 12 ) 

$ = V (®l4 n > m) + <f>2T ( e n > m) ) . (13) 

where the su(2) matrices Tr n ' m \ Tg' m ' , and are defined as products of the spatial unit vectors 

-(«,m) _ (sixi(rn9) cos(tw^), sin(rn6') sxn.(rup), cos(rn6)) , 
^(n,m) _ ( C0S ^ m Q^ CO s(niy9), cos(m#) sm.(nip), — sm(m9)) , 

= (-sin(n(^),cos(n^),0) , (14) 
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with the Pauli matrices r a = (t x ,t v ,t z ), i.e. 

4 n ' m) = 8m(m6)A n) + cos(m6»)r z , 



t^ 1 ' = - sm(rup)r x + cos(rup)T y , 



with Tp = cos(n<p)Tj; +sin(n(p)r 7) . For m = 2, n = 1 the Ansatz corresponds to the one for the monopole-antimonopole 



pair solutions |12lll3Lll8j. while for m = 1, n > 1 it corresponds to the Ansatz for axially symmetric multimonopoles 
Emm The four gauge field functions Ki and two Higgs field functions $i depend on the coordinates r and 9, only. 
All functions are even or odd w.r.t. reflection symmetry, z — > —z. 
The gauge transformation 

U = cxp{iT(r,9)T^/2} (15) 

leaves the Ansatz form-invariant |23j | . To construct regular solutions we have to fix the gauge . Here we impose the 
gauge condition 

rd r K x - d e K 2 = . (16) 

With this Ansatz the equations of motion reduce to a set of 9 coupled partial differential equations, to be solved 
numerically subject to the set of boundary conditions, discussed below. 



C. Boundary conditions 

To obtain globally reg ular asymptotically flat solutions with the proper symmetries, we must impose appropriate 
boundary conditions 0, 0] . 

Boundary conditions at the origin 

Regularity of the solutions at the origin (r = 0) requires for the metric functions the boundary conditions 

dr/M)| r =o = d r m(r,6)\ r=0 = d r l(r,0)\ r =o = , (17) 
whereas the gauge field functions Ki satisfy 

JTi(O,0) = K 3 (O,0) = , K 2 (0,6)=K 4 {0,9) = 1 , (18) 
and the Higgs field functions $^ satisfy 

sin(m(9)$i(0, 6) + cos(m6»)$ 2 (0, 9) = , (19) 



d r [cos(m0)$i(r,0) - sin(m0)$ 2 (r, 9)]\ r=Q = , (20) 

i.e. $,(0,0) =0, <9 r $ z (0,f?) = 0. 
Boundary conditions at infinity 

Asymptotic flatness imposes on the metric functions of the solutions at infinity (r = oo) the boundary conditions 

/ — > 1 , to — > 1 , I — ► 1 . (21) 

Considering the gauge field at infinity, we require that solutions in the vacuum sector Q = 0, where to = 2k, tend to 
a gauge transformed trivial solution, 

e 

and that solutions in the sector with topological charge n, where m = 2k + 1, tend to 
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where 

(n) (l,n) 

$(i,n) = VT (i,n) ^ Aj^W = -^-d9 -nsin6-^—dtp 

is the asymptotic solution of a charge n multimonopole, and U — exxy{—ik9T^}, both for even and odd m. 
In terms of the functions K\ — K4, $1, $2 these boundary conditions read 

K x — ► , K 2 — >l-m, (22) 

cos9 — cos(m#) , , 1 — cos(m6) . . 

A 3 — ► ^ m odd , A 3 — ► ^ — i to even , (23) 

sm fc> sm 9 

AWl-^f, (24) 

$1 — ► 1 , *2 — * . (25) 
Boundary conditions along the symmetry axis 

The boundary conditions along the z-axis (9 — and 9 = tt) are determined by the symmetries. The metric 
functions satisfy along the axis 

dgf = dgm = dgl = , (26) 

whereas the matter field functions satisfy 

Ki = K 3 = $2 = 0, d e K 2 = deKi = d e $i = . (27) 



III. RESULTS AND DISCUSSION 



We have constructed numerically gravitating chains and vortex rings, subject to the above boundary conditions. 
We first briefly address the numerical procedure, and then present our results for the chains (n = 1, 2) and vortex 
rings (n > 3). In particular, we consider the dependence of the solutions on the value of the coupling constant a, and 
restrict to vanishing Higgs self-coupling A = 0. 



A. Numerical procedure 



To construct solutions subject to the above boundary conditions, we map the infinite interval of the variable r onto 
the unit interval of the compactified radial variable x S [0 : 1], 

r 



i.e., the partial derivative with respect to the radial coordinate changes according to 

d r - (1 - X) 2 8 s . 

The numerical calculations are performed with the help of the FIDISOL package based on the Newton-Raphson 
iterative procedure Q • It is therefore essential for the numerical procedure to have a reasonably good initial configura- 
tion. (For details see description and related documentation p3|.) The equations are discretized on a non-equidistant 
grid in x and 9 with typical grids sizes of 70 x 60. The estimates of the relative error for the functions are of the order 
of 10~ 4 , 10~ 3 and 10 -2 for solutions with to = 2, m = 3, 4 and m = 5, 6, respectively. 
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B. Dimensionless quantities 

Let us introduce the dimensionless coordinate x, 

x = r/er . (28) 
The equations then depend only on the dimensionless coupling constant a, 

a 2 = AttGt] 2 . (29) 
The mass M of the solutions is related to the dimensionless mass fj,, via 

"=4^ M < (30) 



where fj, is determined by the derivative of the metric function / at infinity [16l |2 

1 

2^ 



M = _L i im x 2 dxf . (31 ) 



C. Gravitating chains: n — 1 

Let us consider first monopole-antimonopole chains composed of poles of charge n — 1. The topological charge of 
these chains is either zero (for even m) or unity (for odd m). The monopole-antimonopole chains possess m nodes of 
the Higgs field on the z-axis. Because of reflection symmetry, to each node on the positive z-axis there corresponds a 
node on the negative z-axis. For even m (m = 2k) the Higgs field does not have a node at the origin, while for odd m 
(m = 2k + 1) it does. Associating with each pole the location of a single magnetic charge, these chains then possess 
a total of m magnetic monopoles and antimonopoles, located in alternating order on the symmetry axis [l4j . 

The dependence of the m = 1 solution, i.e. of the 't Hooft-Polyakov monopole, on gravity has been studied before. 
When gravity is coupled, a branch of gravitating monopoles emerges from the flat space solution and extends up to 
a maximal value a max , where it merges with a second branch. This second branch extends slightly backwards, and 
bifurcates at a critical value a cr with the branch of extremal RN solutions of unit charge. 

The m = 2 chain, i.e. the monopole-antimonopole pair, shows a different dependence on the coupling constant 
|13|. Again, when gravity is coupled, a branch of gravitating monopole-antimonopole pair solutions emerges from the 
flat space solution, and merges with a second branch of monopole-antimonopole pair solutions at a maximal value of 
the coupling constant a max . But since the monopole-antimonopole pair solutions reside in the topologically trivial 
sector and thus carry no magnetic charge, this upper branch of solutions cannot bifurcate with a branch of extremal 
RN solutions at a critical value a cr of the coupling constant. Instead the upper branch extends all the way back to 
a = 0. Along the first branch the mass of the solutions decreases with increasing a, since with increasing gravitational 
strength the attraction in the system increases. Along the second branch, in contrast, the mass increases strongly 
with decreasing a, and the solutions shrink correspondingly. In the limit of vanishing coupling constant the mass 
then diverges and the solutions shrink to zero size. 

Scaling the coordinates and the Higgs field, however, via 

x = - , $ = a$ , (32) 
a 

leads to a limiting solution with finite size and finite scaled mass fi |13| . 

fi = a/i . (33) 

Indeed, after the scaling the field equations do not depend on a. Instead a appears in the asymptotic boundary 
conditions of the Higgs field, |$| — > a. The Higgs field then becomes trivial on the upper branch as a — ► 0. 
Consequently, the scaled monopole-antimonopole solutions approach an EYM solution as a — > 0. This limiting 
solution is the lowest mass BM solution |l9|. The mass and the scaled mass of the monopole-antimonopole pair 
solutions are exhibited in Fig. 1. 
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(a) 01 (b) a 

Fig. 1 The mass fi (a) and the scaled mass jl (b) are shown as functions of the coupling constant a for the chain solutions with 
n — 1 and m = 2 — 6. The dotted lines extend the EYMH curves of the scaled mass to the mass of the lowest Bartnik-McKinnon 
solution. 

Let us now consider chains with even m, m > 2, which, like the m = 2 monopole-antimonopole pair, also reside in 
the vacuum sector. As one may expect, these chains exhibit an analogous dependence on the coupling constant a as 
the monopole-antimonopole pair. In the limit a — > 0, a lower branch of gravitating solutions emerges from the flat 
space solution, and merges at a maximal value a max with an upper branch of solutions, which extends all the way 
back to a — > 0. The mass of the chains with m = 4 and m = 6 is exhibited in Fig. la. Here, for a — > 0, the emergence 
of the lower branches from the YMH limit is seen as well as the divergence along the upper branches. In Fig. lb the 
scaled mass is exhibited. Clearly, the scaled mass of these solutions approaches the mass of the lowest BM solution in 
the limit a — along the upper branch. The values of the metric functions / and I at the origin are shown in Fig. 2 
for both branches of the chains with m = 4 and m = 6. They evolve from unity in the flat space limit and reach the 
corresponding values of the lowest mass BM solution, when a — ► on the upper branch. 



n=1 n=1 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.1 0.2 0.3 0.4 0.5 0.6 0.7 



(a) 01 (b) « 

Fig. 2 The values of the metric functions / (a) and / (b) at the origin are shown as functions of the coupling constant a 
for the chain solutions with n = 1 and m = 2 — 6. The dotted lines extend the EYMH curves to values of the lowest mass 
Bartnik-McKinnon solution. 

Let us now consider the nodes of the Higgs field for these even m chains. In Fig. 3a we recall the dependence of 
the single positive node of the monopole-antimonopole pair on the coupling constant. Starting from the flat space 
limit, the location of the node continuously moves inwards along both branches, and reaches the origin in the a — > 
limit on the upper branch. (We do not consider excited monopole-antimonopole pair solutions with several nodes 
here, which are related to excited Bartnik-McKinnon solutions |13|.l The flat space m — 4 solution has two positive 
nodes. The evolution of these nodes along both branches is exhibited in Fig. 3b. The location of the inner node again 
moves continuously inwards along both branches, and reaches the origin in the limit a — > on the upper branch. 
In contrast, the location of the outer node reaches a finite limiting value in the limit a — > on the upper branch, 
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which, interestingly, agrees with the location of the single node of the flat space monopole-antimonopole pair solution. 
Similarly, as seen in Fig. 3c, of the 3 positive nodes of the to = 6 solution the location of the innermost node reaches 
the origin in the limit a^Oon the upper branch, whereas the locations of the other two nodes reach finite limiting 
values, which agree with the locations of the two nodes of the flat space to = 4 chain. 



n=1,m=2 n-1.m-4 n=1,m=6 




(d) « (e) « 

Fig. 3 The nodes of the Higgs field are shown as functions of the coupling constant a for the chain solutions with n = \ and 
m = 2 (a), m = 4 (b) m = 6 (c), m = 3 (d), m = 5 (e). The dotted lines extend the nodes of the EYMH solutions to zero or 
to the respective values of the m — 2 YMH solutions. 

To obtain a better understanding of the limit a — > on the upper branch for these chains, let us consider the 
matter and metric functions for small values of the coupling constant. In Fig. 4, for instance, we exhibit the gauge 
field function K2, the modulus of the Higgs field |$| and the metric function / for the to = 4 chain for a small value 
of the coupling constant, a 2 = 0.001, on the upper branch. We observe, that we have to distinguish two regions for 
the matter function K2, an outer region and an inner region. In the outer region, the gauge field function K2 of 
the m = 4 chain agrees well with the (by 2 shifted) corresponding gauge field function K2 of the flat space (to = 2) 
monopole-antimonopole pair solution, while in the inner region the gauge field function K 2 of the to — 4 chain agrees 
well with the corresponding function of the BM solution (in by a scaled coordinates). Likewise, in the outer region 
the modulus of the Higgs field |<f>| of the to = 4 chain agrees well with the modulus of the Higgs field of the flat space 
(to = 2) monopole-antimonopole pair solution. In the BM solution, of course, no Higgs field is present. (An analogous 
pattern holds for the functions K\, K3, K4, $1 and $2 after a gauge transformation.) The metric function /, on the 
other hand, agrees well with the metric function of the BM solution (in by a scaled coordinates) everywhere. (The 
same holds for the functions to and /.) The metric of the flat space (m = 2) monopole-antimonopole pair solution is, 
of course, trivial. 

Thus, for small a on the upper branch, the to = 4 chain may be thought of as composed of a scaled BM solution in 
the inner region and a flat space (to = 2) monopole-antimonopole pair solution in the outer region. Taking the limit 
a — > then, the inner region shrinks to zero size, yielding a solution which is singular at the origin and which has 
infinite mass, while elsewhere it corresponds to the flat space (to = 2) monopole-antimonopole pair solution. This 
explains the above observation concerning the location of the Higgs field nodes, made in Fig. 3. On the other hand, 
taking the limit a — > with scaled coordinates, the inner region covers all of space and contains the BM solution with 
its finite scaled mass, while the fiat space (to = 2) monopole-antimonopole pair solution is shifted all the way out to 
infinity. 
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n=1, m=4 n=1, m=4 




n=1, m=4 




Fig. 4 The gauge field function K2 (a), the modulus of the Higgs field |$| (b) and the metric function / (c) are shown for 
the chain solution with n = 1 and m = 4 at a 2 = 0.001. Also shown are the functions K2 and |$| of the flat space monopole- 
antimonopole solution, and the functions K2 and / of the lowest mass Bartnik-McKinnon solution (in scaled coordinates). 

For the m = 6 chain, we observe an analogous pattern for the functions in the limit a — > on the upper branch, 
as exhibited in Fig. 5. In the outer region, the gauge field function K 2 of the m = 6 chain agrees well with the (by 2 
shifted) corresponding gauge field function Ki of the flat space m — 4 chain, while in the inner region it agrees well 
with the corresponding BM function (in by a scaled coordinates). Likewise, in the outer region the modulus of the 
Higgs field |$| of the m = 6 chain agrees well with the modulus of the Higgs field of the flat space m = 4 chain, while 
the metric function / of the m = 6 chain agrees well with the metric function of the BM solution (in by a scaled 
coordinates) everywhere. 
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Turning now to monopole-antimonopolc chains in the topological sector with unit charge, one might expect a 
completely different pattern for the coupling constant dependence of the solutions, namely a pattern which would be 
similar to the coupling constant dependence of the gravitating monopoles |l5j , where the upper branch of gravitating 
monopoles bifurcates with the branch of extremal RN solutions [HI 113 ■ Such an expectation does not prove true, 
however. Instead, the odd m chains show a similar pattern as the even m chains. Two branches of gravitating m 
chain solutions merge at a maximal value a max . The lower branch of solutions emerges from the flat space m chain, 
and the upper branch extends all the way back to a — 0. For small values of a on the upper branch, the solutions 
may be thought of as composed of a scaled BM solution in the inner region and a flat space m — 2 solution in the 
outer region. 

For the m = 3 and m = 5 chains, we exhibit the mass and scaled mass in Fig. 1, the values of the metric functions 
/ and / at the origin in Fig. 2, and the location of the positive node(s), in Fig. 3. Note, that in the limit a — > on 
the upper branch, the location of the outer positive node of the m = 5 chain agrees with the location of the positive 
node of the flat space m = 3 chain. The composition of the solutions in terms of a scaled BM solution in the inner 
region and a flat space m — 2 solution in the outer region is illustrated in Fig. 6, where we exihibit the gauge field 
function K2 for the m — 3 and m = 5 chains. For the m = 3 chain clearly the gauge field function corresponds in 
the outer region to the function of a flat space (m = 1) monopole, while for the m = 5 chain it corresponds to a flat 
space m = 3 chain. 



n=1, m=3 



n=1, m=5 
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(b) 



Fig. 6 The gauge field function K2 is shown for the chain solutions with n — 1 and m — 3 at a — 0.025 (a) and with n = 1 and 
m — 5 at a 2 — 0.001 (b). Also shown are the functions Ki of the flat space monopole and m = 3 solutions, and the function 
K2 of the lowest mass Bartnik-McKinnon solution (in scaled coordinates). 



Concluding, we observe the same general pattern for all gravitating chain solutions with n = 1. From the flat space 
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to chain a lower branch of gravitating m chains emerges, which merges at a maximal value a max with an upper branch. 
The value of a max decreases with increasing to. For small a the solutions on the upper branch may be thought of as 
composed of a scaled BM solution in the inner region and a flat space to — 2 solution in the outer region. Consequently, 
the mass diverges in the limit a — > 0, while the scaled mass approaches the mass of the lowest BM solution of EYM 
theory. 



D. Gravitating chains: n = 2 

For the m chains composed of monopoles and antimonopoles with charge n = 2 0, [2{| , we observe a completely 
analogous pattern as for the chains composed of monopoles and antimonopoles with charge n — 1. For small a the 
solutions on the upper branch may be thought of as composed of a scaled generalized BM solution with n = 2 in the 
inner region po|. and a flat space m — 2 solution in the outer region. Consequently, the mass again diverges in the 
limit a — * 0, while the scaled mass approaches the mass of the generalized BM solution with n = 2. 

We illustrate this behaviour in Fig. 7 for the m = 4 chain with n — 2. Clearly, the upper branch of the scaled mass 
approaches the mass of the generalized BM solution with n = 2, and the values of the metric functions / and I at 
the origin approach those of the generalized BM solution with n = 2 as well. The composition of the solutions on the 
upper branch for small values of a in terms of an inner generalized BM solution with n = 2 and an outer flat space 
to = 2 solution is exhibited in Fig. 8 for the gauge field function Ki , the modulus of the Higgs field $ and the metric 
function /. 

m=4 m=4 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.1 0.2 0.3 0.4 0.5 0.6 0.7 

(a) « (b) a 

Fig. 7 The mass /i (a) and the scaled mass fi (b) are shown as functions of the coupling constant a for the chain solutions with 
n — 1, 2, 4 and m = 4. The dotted lines extend the EYMH curves of the scaled mass to the masses of the lowest (generalized) 
Bartnik-McKinnon solutions. 
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Fig. 8 The gauge field function K2 (a), the modulus of the Higgs field |$| (b) and the metric function / (c) are shown for the 
chain solution with n = 2 and m = 4 at a = 0.01. Also shown are the functions K2 and |$| of the flat space m — 2 chain, and 
the functions K2 and / of the generalized BM solution with n = 2 (in scaled coordinates) . 



E. Gravitating vortex rings: n = 3 



For larger values of n the flat space solutions completely change character, since the structure of the nodes of the 
Higgs field in these solutions is totally different 0]. Whereas m chains (with n < 2) possess only isolated nodes on 
the symmetry axis, for the solutions with n > 3 the modulus of the Higgs field vanishes on rings in the xy-plane or in 
planes parallel to the a;j/-plane, centered around the z axis. For even m the solutions possess only such vortex rings 
and no nodes on the symmetry axis. For odd m they possess vortex rings as well as a node at the origin, where a 
charge n-monopole is located. 

Let us now consider the effect of gravity on vortex solutions with n = 3. Here the same general pattern is seen for 
the coupling constant dependence as for the gravitating chains. From the flat space vortex solution a lower branch 
of gravitating vortex solutions emerges, and merges at a maximal value a max with an upper branch. On the upper 
branch, for small values of a the solutions may be thought of as composed of a scaled generalized BM solution with 
n = 3 in the inner region and a flat space m — 2 solution in the outer region. Thus the mass diverges in the limit 
a — > on the upper branch, while the scaled mass approaches the mass of the generalized BM solution with n = 3. 
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F. Gravitating vortex rings: n = 4 

It is now tempting to conclude, that for arbitrary values of m and n the same general pattern holds: From the 
flat space YMH solution a lower branch of gravitating EYMH solutions emerges, which merges at a maximal value 
a max with an upper branch. On the upper branch, for small values of a the solutions may be thought of as composed 
of a scaled generalized BM solution in the inner region and a flat space m — 2 solution in the outer region. This 
conjectured pattern is indeed observed for the m — 2 and m = 3 EYMH solutions with n = 4. 

For n = 4 and m = 4, however, a surprise in encountered. As illustrated in Fig. 7, on the upper branch the scaled 
mass of the n = 4 and m = 4 solutions does not approach the mass of the generalized BM solution with n = 4 [2(J , but 
instead it approaches the mass of a new axially symmetric n = 4 EYM solution [2l| . Moreover, a second new axially 
symmetric n = 4 EYM solution exists, which is slightly higher in mass 21] . This second EYM solution constitutes 
the cndpoint of a second upper branch of n = 4 and m = 4 EYMH solutions, which merges at a second maximal value 
of the coupling constant a with a second lower branch, and it is along this second lower branch, that for small values 
of a the solutions may be thought of as composed of a scaled generalized BM solution with n = 4 in the inner region 
and a flat space m — 2 solution in the outer region. Thus for n = 4 and m = 4 there are four branches of solutions 
instead of two, and consequently there are four limiting solutions when a —> 0. The masses of the generalized BM 
solutions and of the new EYM solutions are exhibited in Fig. 9 |2l|. Note, that the boundary conditions of the new 
EYM solutions differ from those of the generalized BM solutions at infinity |2l|. The EYM solutions can thus be 
classified by these boundary conditions. Denoting m = 2k, the new solutions have k = 2, whereas the generalized BM 
solutions have k = 1 |2lj . 




n 



Fig. 9 The dependence of the mass fi on n is shown for the generalized BM solutions (k — 1) and the new EYM solutions 
(k = 2, 3). 

In Fig. 10 we illustrate the dependence of the location of the vortex of the Higgs field on the coupling constant for 
these four branches of n = 4 and m = 4 solutions. The flat space solution has two vortex rings, located symmetrically 
in planes parallel to the xy plane. Starting from the flat space location, the vortex rings move continuously inwards 
and towards the xy plane along the first lower branch and the first upper branch, until at a critical value of the 
coupling constant the two rings merge in the xy plane. When a is decreased further along the first upper branch, a 
bifurcation takes place and two distinct rings appear in the xy plane. In the limit a — > they shrink to zero size. In 
scaled coordinates, in contrast, they reach a finite limiting size. Likewise, in scaled coordinates there are two distinct 
rings with finite size on the second upper branch (located not too far from the two rings of the first upper branch). 
Following the second upper branch and then the second lower branch, the inner ring shrinks to zero size in the limit 
a — > (in ordinary coordinates), while the outer ring approaches the location of the flat space vortex ring of the 
n = 4 and m — 2 solution. 
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Fig. 10 The coordinates p and z (a) and the scaled coordinates p and z (b) are shown as functions of the coupling constant 
a for the m — 4 and n = 4 vortex solutions along the four branches, po, zo refer to the location of the rings of the first two 
branches, and pi , p2 to the location of the rings of the second two branches. 

Summarizing, we find the following pattern: From the fiat space m = 4 and n = 4 vortex solution the first lower 
branch of gravitating vortex solutions emerges, which merges at the first maximal value cn max with the first upper 
branch. Along the first upper branch, for small values of a the solutions approach the scaled k = 2 and n = 4 EYM 
solution with lower mass. Close to the first upper branch, for small values of a the second upper branch arises from 
the scaled k — 2 and n = 4 EYM solution with higher mass. This second upper branch merges at the second smaller 
maximal value a max with the second lower branch of solutions. Along the second lower branch, for small values of a 
the solutions may be thought of as composed of a scaled k = 1 and n — 4 EYM solution in the inner region and a flat 
space m — 2 and n = 4 solution in the outer region. 

As seen in Fig. 9, no EYM solutions with k > 2 and n = 4 appear. We therefore conjecture, that for solutions with 
m > 4 and n = 4 an analogous but generalized pattern holds. There are four branches of solutions. From the flat 
space m > 4 and n — 4 vortex solution the first lower branch emerges, and merges at a first maximal value a max with 
the first upper branch. Along the first upper branch, the solutions approach a composite solution for small values of 
a, consisting of the scaled k = 2 and n = 4 EYM solution with lower mass in the inner region and an m — 4 and 
n = 4 YMH solution in the outer region. The second upper branch arises from a composite solution consisting of the 
scaled k — 2 and n — 4 EYM solution with higher mass in the inner region and an m — 4 and n — 4 YMH solution 
in the outer region. The second upper branch merges at a second maximal value a max with the second lower branch 
of solutions. Along the second lower branch, for small values of a the solutions may be thought of as composed of 
a scaled k = 1 and n = 4 EYM solution in the inner region and a flat space m — 2 and n — 4 solution in the outer 
region. 



G. Gravitating vortex rings: n > 4 



For vortex solutions with n = 5we expect an analogous pattern as for the n — 4 vortex solutions. In particular, 
there should be four branches of solutions: a first lower branch emerging from the flat space solution, a first upper 
branch approaching a composite solution consisting of the lower scaled k = 2 EYM solution and a m — 4 YMH 
solution, a second upper branch arising from a composite solution consisting of the higher scaled k = 2 EYM solution 
and a m — 4 YMH solution, and a second lower branch approaching a composite solution consisting of the scaled 
k = 1 EYM solution and a m — 2 YMH solution. We illustrate this pattern for the m = 5 and n — 5 vortex solutions 
in Fig. 11, where we show the gauge field function K 2 for solutions on all four branches, for small values of a. The 
scaled mass along the four branches is shown in Fig. 11c. 



15 



n=5, m=5 



n=5, m=5 







— e=o 






- - 6=7l/2 


upper branch \ 


y lower branch - 


\ 


- \ 





(a) 



-1 o 

iog 10 M 



1 2 







9=0 






- - 8=71/2 


\ 1 >. 

\ 


-c/ \ upper branch - 


lower branch \ J^r~ 








\ J 





(b) 



-2 -1 

iog 10 (x) 




Fig. 11 The gauge field function K2 is shown for the vortex solution with n — 5 and m = 5 for a 2 — 0.001 on the first lower 
branch and the first upper branch (a), on the second lower branch and the second upper branch (b). The scaled mass /} is 
shown as function of the coupling constant a for all four branches (c). The dotted lines extend the EYMH curves of the scaled 
mass to the masses of the k = 2 and k = 1 EYM solutions. 



Considering higher values of n, we observe in Fig. 9, that for n = 6 and m > 6 at least five EYM solutions exist: one 
k = 1, two k = 2 and two k — 3 solutions. This suggests, that for n — 6 and m = 6 six branches of gravitating vortex 
solutions exist: a first lower branch emerging from the n — 6 flat space solution, a first upper branch approaching the 
lower mass scaled k = 3 and n = 6 EYM solution, a second upper branch arising from the higher mass scaled k = 3 
and n = 6 EYM solution, a second lower branch approaching a composite solution consisting of the lower mass scaled 
k = 2 and n = 6 EYM solution and a m — 4 and n = 6 YMH solution, a third upper branch arising from a composite 
solution consisting of the higher mass scaled k = 2 and n = 6 EYM solution and a m — 4 and n — 6 YMH solution, 
and a third lower branch approaching a composite solution consisting of the scaled k = 1 and n = 6 EYM solution 
and a vn — 2 and n = 6 YMH solution. We indeed obtain these six branches numerically. Surprisingly, there are two 
more branches for a small range of a. We exhibit the scaled mass along these branches in Fig. 12a. The value of the 
metric function I at the origin is exhibited in Fig. 12b. 
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Fig. 12 The scaled mass jl (a) and the value of the metric function I at the origin (b) are shown as functions of the coupling 
constant a for the six branches of the vortex solutions with n — 6 and m = 6. The dotted lines extend the EYMH curves to 
the values of the k — 3, k — 2 and k — 1 EYM solutions. 



IV. CONCLUSIONS 



We have considered the effect of gravity on the monopole-antimonopole chains and vortex rings of Yang-Mills- 
Higgs theory [l4|. The resulting solutions are regular, static, axially symmetric and asymptotically flat. They are 
characterized by two integers m and n. Solutions with n = 1 and n = 2 correspond to gravitating chains of m 
monopoles and antimonopoles, each of charge ±n, where the Higgs field vanishes at m isolated points along the 
symmetry axis. Larger values of n give rise to gravitating vortex solutions, where the Higgs field vanishes on one or 
more rings, centered around the symmetry axis. 

Concerning their dependence on the coupling constant a, we observe the same general pattern for all gravitating 
chain solutions. From the flat space m chain a lower branch of gravitating m chains emerges, which merges at a 
maximal value a max with an upper branch. For small a on the upper branch, the solutions may be thought of as 
composed of a scaled (generalized) BM solution in the inner region and a flat space m — 2 solution in the outer 
region. Consequently, the mass diverges in the limit a —> 0, while the scaled mass approaches the mass of the lowest 
(generalized) BM solution of EYM theory. 

The same pattern holds for gravitating vortex solutions with n = 3, and for n > 4 when m < 3. For n = 4 and 
rn > 4, however, a new phenomenon occurs. Instead of two branches of solutions four branches of solutions arise. The 
reason for these is the presence of two additional EYM solutions [2l|. In this case we find (and conjecture for higher 
m) the following pattern for the coupling constant dependence: A first lower branch emerges from the flat space 
solution, a first upper branch approaches a composite solution consisting of the lower scaled k = 2 EYM solution and 
a m — 4 YMH solution, a second upper branch arises from a composite solution consisting of the higher scaled k = 2 
EYM solution and a m — 4 YMH solution, a second lower branch approaches a composite solution consisting of the 
scaled k = 1 EYM solution and a m — 2 YMH solution. 

For n > 6 and m > 6 not only two but four additional EYM solutions arise [2l| . We consequently find (and 
conjecture for higher m) six branches of gravitating vortex solutions, and for a small range of a even eight branches. 
For larger values of n we expect this trend to continue. 

We observe that for all the gravitating monopole-antimonpole pair solutions, chains of monopoles and antimonopoles 
and vortex rings we considered in this paper there is no formation of a degenerate horizon, in contrast to the gravitating 
monopole solution [l5| . 

We note that for gravitating Skyrmions 29] the dependence of the solutions on the gravitational coupling parameter 
follows the same pattern as for the monopole-antimonpole pair. Again, two branches of solutions merge at the maximal 
value of the coupling parameter, and the upper branch is connected to the (scaled) lowest mass BM solution. 

In this paper we have only considered gravitating chain and vortex solutions in the limit of vanishing Higgs self- 
coupling. For finite values of the self-coupling constant the structure of the flat space solutions can become more 
involved, yielding solutions with vortex rings and several isolated nodes of the Higgs field 14]. Their continuation to 
curved space remains to be studied. 

Furthermore, it would be interesting to consider dyonic gravitating chain and vortex solutions (2(|, as well as 
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black holes with monopole or dipole hair, consisting of chain-like or vortex-like structures [27J. Finally, recent work 
on gravitating monopoles and monopole-antimonopole pairs in the presence of a cosmological constant might be 
generalized to the study of gravitating chain and vortex solutions ■ 
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